Quantum Phase Transition in a Pseudo-hermitian Dicke model 
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We show that a Dicke-type non-hermitian Hamiltonian admits entirely real spectra by mapping it 
to the "Dressed Dicke Model" (DDM) through a similarity transformation. We find a positive-definite 
metric in the Hilbert space of the non-hermitian Hamiltonian so that the time-evolution is unitary 
and allows a consistent quantum description. We then show that this non-hermitian Hamiltonian 
describing non-dissipative quantum processes undergoes Quantum Phase Transition(QPT). The 
exactly solvable limit of the non-hermitian Hamiltonian has also been discussed. 
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Although the choice of a proper set of hcrmitian opera- 
tors is sufficient to ensure the reality of the entire spectra 
and unitary time-evolution for a quantum system, it is 
neither necessary nor dictated by any fundamental prin- 
ciple. It is known since the pioneering work of Bender 
and Boettcher[l| that PT-symmetric non-hermitian op- 
erators with an appropriate inner-product in the Hilbert 
space give consistent description of non-dissipative quan- 
tum processes. The Hamiltonian that is non-hermitian 
with respect to the conventional inner-product in the 
Hilbert space becomes hermitian with respect to the 
new inner product and results of a hermitian theory fol- 
low naturally. The same problem can be studied using 
pseudo-hermitian operator^, 0|, i.e. an operator that 
is related to its adjoint through a similarity transforma- 
tion. Both the approaches involving pseudo-hermiticity 
and "PT-invariance are complementary to each other and 
open up several new directions in the study of non- 
hermitian operators 0, 0, 0, H 0, 0, 0, 0, 0, EH EI H 



13, 14]. It may be mentioned here that operators which 



are non-hermitian with respect to the conventional inner 
product in the Hilbert space are generally used to simu- 
late dissipative processes. In this letter, we are concerned 
about a subclass of such non-hermitian operators which 
are also pseudo-hermitian and may be used consistently 
to describe non-dissipative processes with a modified in- 
ner product in the Hilbert space. 

The study on QPT[15] has received considerable atten- 
tion in recent times and reveals many aspects that are 
qualitatively different from that of phase transition at 
finite temperature. The investigations so far are mainly 
restricted to hermitian Hamiltonian, since an entirely real 
spectra with a well-defined ground-state is not guaran- 
teed a priory for a non-hermitian Hamiltonian. The dy- 
namics of QPT in a closed system is governed by non- 
dissipative terms, since the system at zero temperature 
is already in thermal equilibrium. Unlike the phase- 
transitions at finite temperature, the time-evolution from 
one phase to the other is expected to be unitary for a 



system undergoing QPT. It is thus nontrivial for a non- 
hermitian Hamiltonian to describe a QPT in a closed 
system, since the time-evolution is not necessarily uni- 
tary. It is natural to ask at this juncture whether or not 
one could discuss about QPT in a closed system within 
the framework of PT-symmetric non-hermitian Hamilto- 
nian. If the answer is in the affirmative, it may open up 
new directions in the study of several interlinked areas of 
physics like level statistics, quantum entanglement, quan- 
tum chaos etc. within the framework of pseudo-hermitian 
and/or PT-symmetric non-hermitian Hamiltonian. The 
enlarged parameter-space of a non-hermitian Hamilto- 
nian compared to its hermitian counterpart may prove 
to be an added advantage. 

One of the main results of this letter is that a 
pseudo-hermitian deformation of the DDM indeed un- 
dergoes QPT. We consider the non-hermitian Dicke-type 
Hamiltonian lldl , 



H = waWtfie^V +9 2 e 



J- a 



+ f3e-' 1 ^- J+a + -fe^ 3 J^a + 5e- li3 J+a^ + uj a J z ,{l) 

where aj,u>o,6x, 02, ot,0, 7, <5, ^1,^2,^3 are real parame- 
ters; a, a* are the standard bosonic annihilation-creation 
operators and J z , J± :— J x ± iJ y are the generators of 
the 577(2) algebra, 



[a,a f ] = 1, 
[J+,J_] =2J Z , 



[J Z ,J ± ] =±J±. 



(2) 



The Hamiltonian H commutes with the parity operator 



n, 



n = e' 



riV 



N = a^a + J z +j, 



(3) 



where j is the total spin-angular momentum. The eigen- 
states of H have definite parity depending on whether 
the eigenvalues of the operator N are odd or even. In 
general, the Hamiltonian H is non-hermitian. The Her- 
mitian Hamiltonian is obtained in the limit, 



a = /3, 7 = S, 61 = 



(4) 
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and is known as the DDM in the literature 171 UM ■ The 
standard Dicke model is obtained by a further choice 

of 0i = 2 = & = & = & = and a = P = 
7 = 5. The Dicke Hamiltonian has been studied ex- 
tensively from the viewpoint of QPT 19, 2 
statistics [2 i f. quantum entanglement [22|, |2 
solvability [24j. Certain spintronics based models [24 |25| 
with Dresselhaus and Rashba-type spin-orbit interactions 
can be mapped to the Dicke model, implying its relevance 
in the study of two dimensional semi-conductor physics. 
The Tavis-Cummings model 26] is obtained in the limit 
6>i=6>2 = 7 = 5 = and it reduces to the Jaynes- 
Cummings model[27| if the fundamental representation 
of the SU(2) is used. Non-hermitian versions of both 
the Tavis-Cummings and the Jaynes-Cummings models 
have been studied previously [12]. The Hamiltonian with 
ujo = a = /3 = 7 = <5 = 0is known as the 'Swanson 



model' [ll| in the context of "PT-symmetric quantum me- 



chanics and has been studied in some detail 11, 131. In 



this letter, we study the Hamiltonian H with its full gen- 
erality and show the existence of QPT for certain special 
choices of the parameters. 

The Hamiltonian H can be mapped to a hermitian 
Hamiltonian Ti through a similarity transformation when 
the following relations are satisfied, 

a 6 - P 7 = 0, 0i = 2 = 

a 5 l -(3 7 2 = 0, 0i^O^6i 2 . (5) 

To see this, define an operator p and its inverse as, 
P = e 6 , p-^e- 6 , 

6 = \ln(^a + \ln(g)(J. +j ). (6) 

The operator p is positive-definite and well-defined pro- 
vided the following relations are satisfied, 



7i a t 

r 2 >0 > ^ >0 > h°- 



(7) 



The conditions |^ > and ^ > are sufficient to en- 
sure that p has the desired property. The much more 
stringent condition j7]) is used to make the transformed 
Hamiltonian Ti hermitian. The operator O can be con- 
structed for several special cases as follows: 

6 = i zn GO ata ' ^ >o ' a =p=~f= 5 =^ 

0i = 02 = 0, ^ > 0, 1 > 0; 

p o 

HKI)° to+ M?) (J - +i) ' 

^>0, ^>0, 7 = * = 0; 



O 



\ln{^aU+\ln(T) { J z+ j), 



01 7 

-i>0, j>0, a = = 0. 

02 



(8) 



We will be working within the range of the parameters 
defined by Eq. (JT]) unless mentioned otherwise. Using 
the Baker-Campbell-Hausdorff formula, 



e A Be- A = B+[A, B} + -{A, [A, B]] + -[A, [A, [A, B]]] 



3! 



we find, 

Ti = pHp- 1 



(9) 



uja^a + V0L02 ( e lil a 2 + e^ 1 ^ 



a) 2 ^j + lu J z 



+ Vl$ (e* 3 J-a + J+a f ) , (10) 

when the condition ([5]) is satisfied. Note that Ti is her- 
mitian, since 0102, oc/3 and "fS are positive-definite due to 
the condition ([7]). The Hamiltonian H is quasi-hermitian, 
i.e., related to the hermitian Hamiltonian Ti through a 
similarity transformation. The pseudo-hermiticity of H, 
i.e. = ri + Hri^ r 1 , follows automatically where the met- 
ric rj + is given by r)+ := p 2 . The Hamiltonian H that 
is non-hermitian under the Dirac-hermiticity condition 
becomes hermitian with respect to the modified inner- 
product defined in the Hilbert space as, ({u, v)) v+ '■= 
(u,r)+v). In particular, 

(u\Hv)^(Hu\v), ((u\Hv)), 1+ = ((Hu\v)) v+ . (11) 

Thus, with the modified inner- product, the results of a 
hermitian Hamiltonian follow automatically. 

A comment is in order at this point. The atomic inver- 
sion and the mean photon number are determined by the 
expectation values of the operators J z and a^a, respec- 
tively. Both the operators J z and a'a commute with r] + 
and hence, are hermitian with respect to the modified 



-»£ a 



J 4 



inner- product. However, operators like J X1 J y 



and i(a' — a), which are hermitian with respect to the 
Dirac-hermiticity condition, are no longer hermitian with 
respect to the modified inner product. It may be noted 
here that corresponding to each operator A that is hermi- 
tian with respect to the Dirac-hermiticity condition, the 
operator A := p~ x Ap is hermitian with respect to the 
modified inner product 0]. Consequently, the operator A 
is a physical observable in the Hilbert space of H that 
is endowed with the metric rj + . Following this prescrip- 
tion, a set of SU(2) generators those are hermitian with 
respect to the modified inner-product can be constructed 
as follows: 



Ju . 



J x coshT — iJySinhT 



JyCOshT + iJ x sinhT 



1 fori 
1 = —In — - 

4 \(36 



(12) 
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Similarly, annihilation operator a and its adjoint can 
be obtained as, 



(13) 



The non-hermitian Hamiltonian H can be re-written in 
terms of these operators as, 



H = uja'a 



VflA(e* a 2 + e-^ (a*) 2 ) + uj q J z 



a/3 (V« 2 J_tf + e~^ 2 J+a 
^(e^ 3 J-a + e'^ 3 J+a^ , 



(14) 



where J± := J x ± i J y . 

The hermitian Hamiltonian TL has the form of 
the DDM and has been extensively studied in the 
literature 17 , 18| . In general, the Hamiltonian TL is not 
exactly solvable. Using the Bogoliubov transformation, 



she 



e^sinhO 
coshQ 



(15) 



either the counter-rotating terms J_a, J + a* or the dou- 
ble frequency terms a 2 , a) in the Hamiltonian TL can be 
eliminated with all other terms appearing with renormal- 
ized coupling constants. Both the counter-rotating and 
the double frequency terms can be eliminated simultane- 
ously for fixed values of and cj>, if a constraint involving 
the parameters a, (3, 7, 5, 9\ and 2 is also satisfied. Let 
us choose 4> and as, 



-£1, 9 = tanh 



-if A 



A 



: to - 



4M 2 



* , lo 2 > 40x02, (16) 



so that the double-frequency terms are eliminated from 
TL. It may be mentioned here that the choice of <j> and 9 



= — £1, 9 = tanh 



A 



\ 2\/9\92 J 

A = w+(cj 2 -46» 1 6» 2 )^, uj 2 >A9 1 9 2 , (17) 

also removes the double-frequency terms. However, this 
solution leads to unphysical situations and is discarded 
henceforth. Using the condition {§} and demanding the 
removal of counter-rotating terms, the values of 7, 8 and 
£3 are determined as, 

7 = lt' S = w^ & = * i+ &- ( 18 ) 

The Hamiltonian TL can be expressed in terms of the new 
canonical operators b, w as, 

TL = uj Jz + Cltfb + n + fii (e^ 2 J+b + e l « 2 J_&t) , 



fi 



A 



40!< 



a/3 
20i 0l 



(40i 



40i0 2 - cjA 
•a -wA)5 , 



(19) 



which has the form of the Tavis-Cummings model or the 
DDM in the rotating-wave approximation with the mod- 
ified coupling constants. All these coupling constants 
f2o, and f^i are real and f2 is also positive-definite for 

LO 2 > 401 2 . 

The Hamiltonian 7i in Eq. (fl"9)) is exactly solvable [20]. 
It can be decomposed in terms of two mutually commut- 
ing operators K and L as follows, 



TL = Q,K + Q, x L + n Q , 
K = tfb + J z , 



L 



J + b + e ii2 J_6 f 



w 



Qi 



-J, 



(20) 



The operator K is diagonal for a fixed spin j with the 
eigenvalues of J z as (m — j), m = 0,1, ...2j and that 
of the bosonic number operator wb as n, n = 0, 1, 2 . . .. 
The operator _L and hence, the operator TL can be diago- 
nalized in the basis spanned by the eigenstates of K. Let 
\ n , m ij)H be a complete set of orthonormal eigenstates 
of TL with the eigenvalues E n ^ n .j. The orthonormality 
of \n,m;j)fi is based on the standard inner-product in 
the Hilbert space. The eigenstates of H with the same 
eigenvalues E niTn; j are determined as, 



\n,m;j) H = p 1 \n,m;j) n , 



(21) 



which form a complete set of orthonormal eigenstates 
under the modified inner-product defined in the Hilbert 
space of H. Consequently, the non-hermitian Hamilto- 
nian H is also exactly solvable and admits consistent 
quantum description. 

The expectation value of an operator X in the Hilbert 
space of H is determined as, 



v+ 



(n,m;j\pXp 1 \n,m;j) H - 



(22) 



Both J z and a^a are hermitian with respect to the Dirac- 
hermiticity condition as well as with respect to the mod- 
ified inner product. In particular, both J z and a'a com- 
mute with p, leading to the results: 



((Jz))r,+ = ( n , m ;j\Jz\n,m;j) n , 
((a^a))r, + = (n,m;j\a' i a\n,m;j)H- 



(23) 



Thus, both ((J z )) v+ and ((a'a))^ + are real. However, 
in general, ((J x )) v+ , ((J v ))r,+ , ((a + al)) v+ and ((i(a t - 
a))) n+ are complex, 

({Jx)) v + = coshT(J x ) H + isinhT(J y )n 
{(Jy))v+ = -isinhV{J x ) n + coshT(J y } H 
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«« + « t )) t;+ = (7T i ) ' <«}«+(Y) (^)n 



7 = ±a, 5 = ±/3, 



((at-a)), + =( J) ( a t) w - -M (a) H) (24) 



where (J x )n and {J y )n are real, while (a)« and (a^)-n 
are complex. As discussed before, corresponding to 
each operator A in the Hilbert space of Ti, the phys- 
ical observable in the Hilbert space of H that is en- 
dowed with the metric 77+ is A = p~ l Ap. The expec- 
tation values of these capped operators are real, since 
((A)) n+ — (n,m; j\A\n,m; j)u- Thus, a complete and 
consistent description of the pseudo-hermitian H is al- 
lowed with the proper identification of the physical ob- 
servables, 

The Hamiltonian Ti in Eq. (JX^J) is known to exhibit 
qpt[H I21I ]. Although Ti is exactly solvable, it becomes 
tedious to calculate the eigenspectra for large j. The 
Holstein-Primakoff representation of the SU (2) genera- 
tors 

J = (2j - C f C) ' C, J+ = C f (2j - C f C) ' , J* = C f C-i, 

(25) 

where £, are the bosonic annihilation and creation op- 
erators satisfying [CC^] = 1) ca n be used to study the 
thermodynamic limit j — > 00. It is important to note 
that among the SU (2) generators, only the combination 
J z + j appears in the expressions of the parity operator 
II and the similarity operator p. Consequently, p and 
II are well-defined in the thermodynamic limit j_ 



Following the standard method described in[20j, |2l|, the 



normal phase of Ti can be found to be described in the 
range Ai < = y/QujQ, while the 'super-radiant phase' 
is described in the range Ai > Af, where Ai = y/2jH,i. 
The Hamiltonian H and Ti have the same eigenspectra, 
since they are related to each other through a similarity 
transformation. Moreover, note that the operators O, 
p and p _1 are well-defined in the thermodynamic limit 
j — ► 00. Thus, the Hamiltonian H also undergoes QPT 
with the normal phase described in the range Ai < A^, 
while the 'super-radiant phase' is described in the range 
Ai > Af . The values of the mean photon number and 
the atomic inversion above the critical value \\ can be 
determined as follows: 



r l ((a^a)) 



1 



ft? 



(26) 



This is one of the main results of this letter. 

The hcrmitian Hamiltonian Ti in Eq. (fT4"|) with 6\ = 
2 = £1 = £2 = £3 = and = = ;^§j reduces 

to the 'standard Dicke mo del' which is known to un- 
dergo QPT for A 2 > X c 2 = ;21]. The non-hermitian 
Hamiltonian if in (JTJ with X = 9 2 = 0, £1 = £ 2 = £3 = 0, 



,t 



H = ui a [ a + +ujq J z + a J_a T + (3 J+a±aJ-ai f3J+a 

(27) 

is equivalent to the 'standard Dicke Model' through the 
similarity transformation H Dicke = pHp^ 1 with the op- 
erator O given by, 

= \ln(^)(J z +j), ^>0. (28) 

Thus, the non-hermitian Hamiltonian H also undergoes 
QPT for A 2 > A§. The values of the atomic inversion and 
the mean photon number above the critical value A 2 are 
identical to that of the standard Dicke model: 



r l ((Jz))n + = -[f 



r i((at a )), + = 



A 2 > X c 2 .(29) 



The results for finite j, as quoted in Ref. [2l| for H Dicke, 
are equally applicable for H, since ((J z )) n+ — (Jz)H Dicke 
and ((a)a)) v+ = (a^a) Hmcke . 

The Hamiltonian Ti in Eq. (|14p with its full generality 
also undergoes QPT for | p \ < 1, 



uj (u + 
(A3 + A4) 2 



A 3 = 




A 4 = 



(30) 



Consequently, H with the parameters satisfying the rela- 
tions in Eq. ([5]) also undergoes quantum phase transition 
for j ji I < 1. The values of mean photon number and the 
atomic inversion for p < 1 can be determined as follows: 



r x «aM) 



1+ 



j 1 ((Jz)) V+ =-p, V < 1- 



A3 + A4 

• + 2VW2 



(31) 



The mean photon number vanishes identically and the 
atomic inversion is equal to —1 for p > 1. The QPT in 
the Tavis-Cummings model and the Dicke model appear 
as special cases of the general result described by Eqs. 
@Q) and (ED- 

We have shown that a non-hermitian version of the 
DDM undergoes QPT. This is the first time in the liter- 
ature that QPT for pseudo-hermitian operators has been 
described and definitely broadens the scope of study- 
ing QPT in various other non-hermitian models. For 
the particular case of the pseudo-hermitian DDM, it 
is to be seen whether or not the QPT is related to a 
change in level-statistics and/or cross-over from entan- 
gled to disentangled states, as is the case for the stan- 
dard Dicke Hamiltonian [U Finally, as mentioned 
earlier, the DDM can be mapped to certain spintronics- 
based models 0, 25 [. Our results on QPT can be di- 



rectly extended to such models and may prove to be 
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the testing ground of pseudo-hermitian quantum me- 
chanics through appropriate quantum engineering of two- 
dimensional semiconductor devices. 
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